Glassy dynamics of the inhomogeneous metallic phase in Lai.^Ca^MnOs 
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The anelastic spectra (elastic energy loss and dynamic modulus versus temperature) of 
Lai-^Ca^MnOs have been measured for 0.33 < x < 0.5. A peak in the imaginary part of the 
elastic susceptibility is found slightly below the Curie temperature, whose temperature and fre- 
quency dependences are typical of the dielectric and magnetic susceptibilities of relaxor ferroelectrics, 
spin glasses and other inhomogeneous frustrated systems. The fluctuations can be associated with 
the magnetic clusters which are observed by several techniques or with charge-ordered insulating 
domains in the low temperature inhomogeneous phase. The temperature dependent distribution 
function of the fluctuation energy barriers has been extracted by the analysis of the experimental 
curves. 
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The perovskite manganese oxides Ri-^A^MnOa (R = 
La, Pr, or Nd and A = Ca, Ba, Sr, or Pb) are the 
object of intense study, since the strong coupling be- 
tween their electronic, magnetic, orbital and structural 
degrees of freedom gives rise to a rich .phase diagram 
as a function of doping and temperature!!! In particular, 
Lai_2;Caa;Mn03 in the doping range 0.20 < x < 0.50 un- 
dergoes a phase transition from paramagnetic insulating 
(PI) to metallic ferromagnetic (FM) upon cooling be- 
low room temperature, and exhibit the so called nega- 
tive colossal magnetoresistance (CMR). The original in- 
terpretation of the magneto-conductive phase transition 
by the classical double exchange (DE) mechanism has 
been recently demonstrated to be not quantitatively suf- 
ficient. Therefore a strong charge-phonon polaronic cou- 
pling, associated with the Jahn- Teller distortion of the 
Mn 3+ (0 2 ~)6 octahedra, has been supposedo, which still 
does not supply an exhaustive interpretation of CMR. In 
particular, a strong magnetic correlation between charge 
carriers and local moments causes the self-trapping of 
the carriers in magnetic polaronstj over regions of about 
10-20 A. When a magnetic field is applied, these ferro- 
magnetic clusters overlap leading to carrier derealization 
and consequently to CMR.tffl Such spin polarons have 
been first revealed in non Jahn- Teller systems presenting 
CMR, like dilute magnetic semiconductors, rare earth 
chalcogenides, Mn pyrochlores, and now also in mangan- 
ites around thepCurie temperature Tc by inelastic neu- 
tron scattexing,B small angle neutron scatteringa (SANS) 
and NMR.Q These magnetic polarons can be the precur- 
sors of the phase-.inhomogeneities recently evidenced be^ 
low T c by STM,B NMR0 and X-ray diffraction studies.EHl 
In this phase, insulating and metallic domains of 10- 
50 nm coexist and percolate as a function of tempera- 
ture and applied magnetic field. Such a phase separated 
scenario recalls the relaxor ferroelectrics, which also are 
perovskites, where the mixed composition induces fer- 
roelectric clusters of nanometer dimensions embedded 
in a paraelectric matrixOEij Analogously, in mangan- 
ites with mixed compositions lile Lai^^-yPryCa^MnOa 
and Ndo.sCao.sMni-yCr^Oa it is possible to induce the 
coexistence of domains of FM and charge-orbital or- 
dered (CO-OO) phases, which are stable and visible by 
electron diffraction and by electron and magnetic force 
microscopyJ13o In these cases aging effects in magnetiza- 
tion and conductance have been observed just like those 
revealed in glassy magnetic or relaxor ferroelectrics. It 
seems that in the CMR manganites this phase separation 
is more overbalanced towards the ferromagnetic compo- 
nent and the insulating CO component is more fluctuat- 
ing and less coherent, resulting therefore not visible by 
diffraction. Increasing the Ca doping above 0.50, the CO 
phase becomes long ranged.lil 

In this general context it is of great interest to study 
the dynamics of these inhomogeneous fluctuating do- 
mains. This study is possible by anelastic spectroscopy 
measurements, since the elastic constants and the elas- 
tic energy dissipation coefficient are strictly connected 



to the dynamic response of the system,ll3 and also the 
magnetic fluctuations can affect the anelastic response 
through magnetoelastic coupling. 

The Lai-^Ca^MnOa samples, with x ranging from 
0.33 to 0.5, were prepared from La203, CaC03 and 
Mn02 powders and sintered at 1300-1400 °C several 
times, until the X-ray diffraction spectra indicated that 
all the starting powders had reacted; the ingots were cut 
into bars approximately 45 x 4 x 0.5 mm 3 . The com- 
plex dynamic Young's modulus E (oj, T) = E' + iE" was 
measured as a function of temperature by electrostati- 
cally exciting the 1st and 5th free flexural modes, whose 
frequencies uj/2-k are in the ratios 1 : 13.2 and detect- 
ing the vibration amplitude with a frequency modula- 
tion technique. The relative change of the real part E' 
was measured from the temperature dependence of uj 2 = 
E'/p, where the mass density p varies much less than 
E; the imaginary part is related to the elastic energy 
loss coefficient (or reciprocal of the mechanical Q) by 
Q- 1 (uj, T) = E"/E' = X "/x', where the X = 1/E is the 
elastic susceptibility (compliance), and it was measured 
from the decay of the free oscillations or from the width 
of the resonance peak. 

The anelastic spectrum of Lai-aCa^MnOs with 0.33 < 
x < 0.5 is dominated by the trigonal/orthorhombic 
structural transition at T t and the PI/FM transition 
at Tq. Figure 1 presents the whole anelastic spectrum 
of a sample with x = 0.33, with T t = 721 K and 
Tc = 250 K. The sharpness of the structural transition at 
T t = 721 K is an indication of the good sample homogene- 
ity, while the huge stiffening below Tc = 250 K demon- 
strates the strong interplay between orbital/magnetic 
and lattice phenomena, already observed by ultrasonic 
measurements I 16 !! 1 ! Similar spectra were obtained with 
Ca doping up to 50%. 

In the present paper we are concerned with the region 
of the spectrum just below the PI/FM transition, which 
did't present any par.tirrujar feature in the previous ul- 
trasonic experiments.ESO This part of the spectrum is 
shown in Fig. 2 for the sample with x = 0.33. In analyz- 
ing acoustic data, it is always assumed that the absorp- 
tion can be identified with the imaginary susceptibility; 
in the present case, however, we have to analyze the spec- 
trum just below Tc, where the large and rapid stiffening 
of the modulus has to be taken into account. There- 
fore, instead of Q _1 we plot x"/Xo — Q~ 1 x'/xo, where 
X>,T) = xoKO) MT)] 2 and X o = x'(0,0) is the 
reciprocal of the static Young's modulus at T = 0. Be- 
side a sharp attenuation peak centred near Tc = 250 K, 
also determined from the onset of the stiffening (Fig. 1) 
and from the peak in the electrical resistivity, a broad 
shoulder appears at lower temperature. Such a peak is 
the focus of the following discussion, and was observed, 
although with different amplitudes, in all our samples 
with five Ca doping levels included between 0.33 and 0.5; 
we chose to analyze the two cases in which the back- 
ground subtraction was easier. First of all the absorp- 
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tion connected with the PI/FM transition has to be sub- 
tracted. The XLCak cannot be identified with the critical 
absorption Jilo Q" 1 ocw(T- T c )~ e , which was foundEI 
in Lao.67Cao.33Mn03 at 10 MHz with 9 = 1.55; indeed, 
the amplitude of that peak extrapolated to our frequen- 
cies becomes smaller than 1 x 10 -5 , while the height of 
the peak in Fig. 2 is 5 — 8 x 10 -4 (but is smaller in the 
samples with higher Ca doping). Whatever the precise 
origin of the absorpion around Tc, it can be well fit- 
ted with a Lorentzian centred near Tc, as shown by the 
continuous lines in Fig. 2a. After subtraction of the ab- 
sorption around Tc and of a small linear background, the 
two curves in Fig. 2b are obtained. The resulting peak 
shifts to higher temperature with increasing frequency, 
like a Debye relaxation process with a thermally acti- 
vated relaxation time r, usually following the Arrhenius 
law, r = 2p exp (E/T). The response of such a process 
would bdlj 



X (w, T) ocT 1 (1 - iljt) 1 



(1) 



whose imaginary part is peaked at a temperature T„ 
satisfying the condition 



WT (T m ) = 1 



(2) 



However, in that case the whole dispersion would shift to 
higher temperature with increasing frequency, whereas 
the low temperature side of the peak in Fig. 2b, is in- 
dependent of frequency. In addition, the height of the 
peak increases considerably with temperature, instead 
of decreasing as 1/T. All these features are typical of 
the dispersion of the imaginary susceptibility of soin- 
glasses, mixed ferroelecteic/antiferroelectric systemsa or 
relaxor ferroelectricsJiit3 and therefore the absorption 
peak seems to originate from the glassy dynamics of some 
inhomogeneous phases, for example from the formation 
and disappearance or reorientiation of the minority in- 
sulating nanodomains which are believed to coexist with 
the FM phase below Tc. In order to extract a more 
quantititative signature of the glassy dynamics, we per- 
formed on the present data the type of analysis adopted 
for spin-glasses and their analogues, including the relaxor 
ferroelectrics. 

At first we tested whether the data at different fre- 
quencies can be collapsed on a single curve by means of 
an appropriate scaling variable E(lj,T) related to the 
law of the slowing down of the fluctuations. The idea is 
that the susceptibility can be expressed in terms of su- 
perposition of elementary processes, with a very broad 
distribution g (t, T) of relaxation times r, 

X (w, T) = X (0, T) / d (In r) g (t, T) — L_ (3) 

Jr 1 - «WT 

so that each contributing peak in the imaginary part 
can be approximated as a 5 function, lot / 1 + [lot] 



f 5 (ljt - 1) . Then X " (w, T) ~ fx (0, T) g T) is di- 
rectly related to g through a scaling relationship often 
derived from the Vogel-Fulcher law 



T = T exp[E/ (T-To)] 



(4) 



and the condition (g). Following CourtensJIS we tested 
the data against the expression 



X " = f(T)R(E c -E) 



(5) 



where R (x) = [1 + tanh (bx)] /2 produces the high tem- 
perature dispersion, E c is a constant and E is related 
to lj and T through Eqs. (§) and (|). The func- 
tion / (T) is determined by fitting the frequency inde- 
pendent region of X " , and is found to be well approx- 
imated by / (T) = Aexp(W/T) with W = 4460 K. 
The dark continuous curves of Fig. 2b, obtained with 
t = 5.5 x 10~ 12 s, T = 193 K, E c = 620 K and 
b = 0.011, demonstrate the glassy nature of the relax- 
ation process just below Tc, and the parameter To could 
be interpreted as a freezing temperature; the use of a 
simple Arrhenius law, putting To = 0, yielded a definitely 
worse fit. 

The same analysis has been carried out for a sample 
with x = 0.50 nominally, but with actual Ca doping 
slightly below that value, since only the PI/FM transi- 
tion was detected, and the resistivity curve did not show 
the expected insulating behavior at low temperature. At 
this doping the intensity of the glassy relaxation peak 
was higher, but the background subtraction was less ob- 
vious, due to the presence of another peak at lower tem- 
perature; for this reason, the simultaneous optimization 
of the background and peak parameters was preferred. 
Figure 3 shows the resulting fit with the same criteria 
used for the x = 0.33 case. The resulting parameters are 
W = 1920 K, t = 6.4 x 10~ 8 s, T = 180 K, E c = 272 K 
and b = 0.0124. 

While the above fits demonstrate the close similar- 
ity between the dynamics of the inhomogeneities in 
Lai_ 2; Ca a ;Mn03 below Tc and a freezing process, most 
of the parameters have little physical meaning. In or- 
der to obtain a more physical picture, we also car- 
ried out another type of analysis, which was recently 
adopted to describe, the dielectric susceptibility of re- 
laxor ferroelectrics Ji3 The relaxation process is described 
in terms of a distribution of activation energies g (E, T) — 
Eq 1 exp \{E C — E) /Eq] with lower cutoff at E c and a 
temperature dependent width E (T). The integration of 
the elementary relaxation function Eq. ([!]) over g (E, T), 
assuming the Arrhenius. law r = To exp (E/T), yields a 
very simple expressior£3 of x" , which we generalized in- 
cluding the Vogel-Fulcher law, Eq. M). The result is 



X" - Xo 



with n 



2 cos (tvk/2) 
T-T 



£0 



exp (E c /E ) {coT ) n (6) 
(7) 
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The simple result (0) is valid for n < 1, and therefore 
we put an upper cutoff to n, which however does not af- 
fect the fit for the case x = 0.33, since n approaches 1 
only in the high temperature region where the curves are 
dropping to zero. The temperature dependence of Eq can 
be extracted from the frequency independent part of the 
curves at low temperature, where «Cl and, initially as- 
suming T = for simplicity, x" (T) ^ (tt/2) xoT/E ; as 
already mentioned, in that region x" (T) — Aexp (W/T), 
and therefore we adopted E (T) = cTexp (W/T), where 
c is a constant. Since W = 2000 — 4500 K, we are in the 
regime were the width Eq (T) of the distribution func- 
tion diverges with lowering temperature, as required by 
a freezing process. The gray curves of Fig. 2b are ob- 
tained with W = 4480 K, c = 2.5 x 10~ 9 , E c = 880 K, 
r = 4 x 10~ 14 s and T = 186 K. The fit is satisfactory, 
considering the simplificative hypotheses made for the 
distribution function of the activation energies, g (E, Eq) 
with Eq (T) , and characterizes the spectrum of the fluctu- 
ations. The relaxation rates are close to Tq 1 exp (—E c /T) 
at high temperature, but the spectrum quickly broadens 
with decreasing temperature, consistently with a freezing 
temperature To ~ 186 K. In a sense, the broadening of 
g (E, T) itself should describe the freezing process, with- 
out the need of introducing the Vogel-Fulcher law, but 
the limitation of using a tractable expression for g (E, T) 
makes it useful the introduction of To > for improving 
the fit. For x — 0.5 (Fig. 3) the condition n < 1 for the 
validity of Eq. (^) was shifted to slightly lower tempera- 
ture, making the fit less meaningful for the data at higher 
frequency. Anyway, the case x = 0.5 can be described by 
W > 1900 K, E c > 980 K, r < 2 x 10~ 10 s; the use of 
To > does not improve significantly the fit. 

Let us consider now the nature of the fluctua- 
tions which are observed in the anelastic spectra of 
Lai_2:Ca x Mn03. A first possibility are the FM clus- 
ters which are observed on different scales, with the 
anelastic response driven by strong magnetoelastic cou- 
pling. Indeed, around Tq SANS measurementsQ reveal 
the presence of "magnetic polarons" with sizes of the 
order of tens of A, while in neutron spectroscopy ex- 
periments at x — 0.33 a quasielastic components ap- 
pears, which indicate a diffusive component of the spin 
excitations, coexisting with the usual spin waves of a 
typical ferromagnetE It should be noted that, although 
the present measurements probe the very low frequency 
part of the fluctuation spectrum, the resulting distribu- 
tion of activation energies g (E, T) acquires weight near 
E c ~ 880 K with increasing temperature above 220 K; 
this corresponds to very fast relaxation rates, in the range 
of 10~ 12 s _1 , which are the rates probed by neutron spec- 
troscopy. In addition, the relaxation could be connected 
with the fluctuation of the magnetization in the nano- 
metric FM domains, observed by STMS and suggested 
by PDF XRDtiJ and NMFB experiments. Another possi- 
bility is that we observe the CO insulating clusters which 
coexist with the FM phase, and coyJd consist of diagonal 
Mn 4+ /Mn 3+ stripes, as observerjila in the CO coherent 



phase at x > 0.50. In this case the anelastic response 
may result both from the formation and disappearance 
of the CO state over few lattice constants or from the 
reorientation of the Jahn- Teller distorted stripes. 

Regarding the origin of the glassy dynamics, in some 
cases a broad distribution of the activation energies is 
supposed to-Eeflect a broad distribution of the sizes of 
the clusters.E3 We rather suggest that the main cause 
of the glassy dynamics is the competition of interac- 
tions with different length scales, which may produce 
spatial inhomogcncitips. with well defined patterns but 
with glassy dynamics In fact, the relevant interactions 
in the manganites range from the magnetic exchange to 
the Coulomb and elastic ones. The broadening of the dis- 
tribution function g (r, T) with lowering T can then be 
due to the frustrated nature of the interactions between 
the coexisting insulating and metallic regions. 

In conclusion, the low frequency dynamic susceptibil- 
ity of Lai_ x Ca^Mn03, probed by anelastic spectroscopy 
for 0.33 < x < 0.5, contains the signature of glassy dy- 
namics in the low temperature inhomogeneous ferromag- 
netic phase. A peak in the elastic energy loss is found 
slightly below the Curie temperature, whose temperature 
and frequency dependences are typical of the dynamic 
susceptibilities of other perovskite inhomogeneous frus- 
trated systems like the relaxor ferroelectrics, and of spin 
glasses. The data have been analyzed in order to ex- 
tract the temperature dependent distribution function of 
the energy barriers of the fluctuations characterizing the 
inhomogeneous state. 
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FIG. 1. Elastic energy loss coefficient, E"/E', and nor- 
malized real part of the dynamic Young's modulus of 
Lao.67Cao.33MnC>3. 

FIG. 2. Imaginary part of the elastic susceptibility of 
Lao.67Cao.33Mn03 measured at two vibration frequencies. 
The continuous lines in (a) have been subtracted to obtain 
the peak in (b) . The black curves in (b) are the fit according 
to Eq. (^), while the gray ones according to Eq. (|). 

FIG. 3. Imaginary part of the elastic susceptibility of 
Lai-zCa^MnOs with x slightly below 0.50. The fit is made 
as described in the text, according to Eq. (^). 
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